Quantum dynamics in high codimension tilings: from quasiperiodicity to disorder 
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We analyze the spreading of wavepackets in two-dimensional quasiperiodic and random tilings as 
a function of their codimension, i.e. of their topological complexity. In the quasiperiodic case, we 
show that the diffusion exponent that characterizes the propagation decreases when the codimension 
increases and goes to 1/2 in the high codimension limit. By constrast, the exponent for the random 
tilings is independent of their codimension and also equals 1/2. This shows that, in high codimension, 
the quasiperiodicity is irrelevant and that the topological disorder leads in every case, to a diffusive 
regime, at least in the time scale investigated here. 
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It is now well established that quasiperiodic order has a 
strong influence on the quantum dynamics of wavepack- 
ets. Indeed, the nature of the eigenstates in quasiperiodic 
systems, which are neither spatially extended (as in pe- 
riodic systems) nor localized (as in disordered systems) 
but critical, is often responsible for a sub-ballistic mo- 
tion. Although most of the studies about this anomalous 
diffusion concern one-dimensional systems such as the Fi- 
bonacci or the Harper chain, there has also been a great 
interest for the, more physical, higher-dimensional ones 
ELQIEQ>I3' However, the parameters that determine 
the characteristics of the long time dynamics, such as the 
diffusion exponent (3, remains misunderstood. 

In this paper, we investigate the quantum dynamics 
of wavepackets in two-dimensional quasiperiodic tilings 
built with De Bruijn grid method ja0|- This construc- 
tion allows us to change easily the codimension of the 
structures that determines the number of different tiles 
and thus, in a sense, its complexity. The codimension 
reveals itself as a fundamental quantity for the dynam- 
ics since the diffusion exponent decreases when the codi- 
mension increases and converges towards 1/2 in the high 
codimension limit. We also analyze the influence of pha- 
son flips that consists in turning the quasiperiodic tiling 
into a random one without changing the stcechiometry of 
the tiles (see Fig. EJ. We show that contrary to what was 
claimed in Ref. jl| for the octagonal tiling, the phason 
disorder do not fasten the spreading but, always slow it 
down towards a diffusive regime. 

The De Bruijn grid algorithm 

HQ 

relies on two steps. 

First, we define D families of iVj equally spaced paral- 
lel lines, rotated one from the other by an angle 2tt/D. 
These D x Ni lines form a grid and separate different 
connected regions of the plane called cells of the decom- 
position. Second, we dualize this grid by associating a 
tiling vertex to each cell and connecting any vertices as- 
sociated to two cells separated by a line. 

This algorithm is simpler to implement in two di- 
mensions than the standard Cut-and-Project 0, 




FIG. 1: Left: Quasiperiodic tiling built with the grid method 
with D = 4. Right: Random tiling obtained from the 
quasiperiodic tiling by elementary flips (center). Such trans- 
formations conserve the stcechiometry of the tiles but destroy 
quasiperiodicity and change the coordination number of the 
sites. 



Ilfjf method where tilings are generated from high- 
dimensional spaces. In this latter approach, the codi- 
mension of the tiling is defined by the difference between 
the dimension D of the initial lattice in which the points 
are selected, and the dimension d of the tiling. For two- 
dimensional tilings, the Cut-and-Project method requires 
the analyzis of a (D — 2)-dimensional (perpendicular) 
space whose complexity increases very quickly with D 
so that, in practice, this method is not tractable beyond 
D = 5. By contrast, in the grid method, a tiling of codi- 
mension (D—2) is straightforwardly obtained by choosing 
D families of lines. Note that for D = 2, the grid method 
builds the usual square lattice, and, for D — 3, it builds 
the well-known (periodic) dice lattice that displays inter- 
esting properties when embedded in a magnetic field |lll | . 
Nonperiodic (quasiperiodic) tilings are thus obtained for 
D > 4. The number of lines AT; determines the size of the 
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tiling since the number of sites is N s ~ NfD{D — l)/2. 
However, with this method, different tilings can be ob- 
tained depending on the relative position of each line 
families. Here, we have focussed on a generic family of 
tilings and checked that our results were weakly sensitive 
to this dephasing choice. In addition, to avoid spurious 
effects due to fixed boundary conditions (Ni is fi- 
nite), we have isolated a circular central cluster of about 
10 6 sites, which is quasiperiodic. 

The dynamics is given by a standard tight-binding 
Hamiltonian: 

where \i) denotes a localized orbital on site i and where 
the hopping term tij equals to I if i and j are nearest 
neighbors and otherwise. Of course, the time evolu- 
tion of any wavepacket is directly obtained by diagonal- 
izing H, but unfortunately, exact diagonalizations are re- 
stricted to rather small system size (~ 10000 sites). Since 
we are interested in the long time behavior of the dynam- 
ics that requires large systems, we have used an approx- 
imate method, the Second Order Differencing Scheme 
[l3j |. which consists, for a given initial state |^(0)), in 
writing: 

\tp(t + At)) = \i/)(t - At)) - 2iAtH\ip(t)). (2) 

Practically, a time step At — 0.05 is sufficient to get 
a good accuracy and all our results have been obtained 
with this value. Note that despite the very low order of 
this development for the evolution operator e lHt , this al- 
gorithm is actually an very efficient tool (see for example 

El). 

The time evolution of a wavepacket can be character- 
ized by several observables. Here, we focus on the mean 
square spreading defined by: 

AR*(t) = m)\R 2 \m) [{m\R\m)T, (3) 

where R is the position operator, the origin being taken 
at the center of the cluster. The diffusion exponent (3 is 
defined by the long time behavior of AR 2 (t ) ~ t 2 @ . A pri- 
ori, nothing prevents from getting a different asymptotic 
regime but, actually, AR 2 often seems to be rather well 
described by a power-law at large time. The case (3=1 
obtained for periodic potential defines a ballistic propa- 
gation whereas < j3 < 1, observed in many quasiperi- 
odic systems, corresponds to a sub-ballistic spreading. 
Note that the diffusive motion ((3 = 1/2) is obtained 
in several situations among which the three-dimensional 
Anderson model in the metallic part of the spectrum. A 
crucial issue is to understand the parameters which de- 
termine the value of this exponent. 

As already noticed in the octagonal tiling Q for 
wavepackets initially localized on a single site, (3 strongly 



depends on this site. In the Labyrinth tiling 3] and in 
the generalized Rauzy tilings |^ where energy-filtered 
wavepackets have been studied, (3 has also been shown 
to be energy-dependent. In the present work, we do not 
claim to give a precise characterization of the diffusion 
exponent for each tiling but we wish to investigate the be- 
havior of (3 with respect to the codimension. Therefore, 
we have chosen to consider random phase states initially 
spread over a finite portion of the cluster (typically 1%). 
This leads to an exponent which is "averaged" over the 
total density of states. We have considered several con- 
figurations of the phases and several initial radius, and 
checked that although (3 is different from one state to an- 
other, it fluctuates of less than 5%. In addition, since we 
have finite size systems with open boundary conditions, 
the propagation must be stopped when the wavepacket 
reach the boundaries. As a criterion for this maximum 
time i max . , we have chosen the time for which the presence 
probability integrated over all the border sites reaches 
1/N S . Typically, this leads to i max . ~ 300. As usual, we 
can never ensure that the real asymptotic regime (if it 
exists !) is reached but this time scale certainly gives 
some hints concerning the long time behavior. 

The mean square spreading AR 2 (t) is displayed in Fig. 
El for quasiperiodic tilings with D £ [3,6]. The dice lat- 
tice D = 3 is given as a reference for a ballistic motion 
(J3 = 1). For low codimension tilings, we observe a su- 
perdiffusive regime with an exponent that decreases when 
the codimension increases. 
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FIG. 2: Log- log plot of AR 2 (t) for low codimension quasiperi- 
odic tilings. 

More precisely, we get (3 = 0.77(1) for the octagonal 
tiling [D =A),(3 = 0.71(1) for the Penrose tiling [D = 5), 
and (3 = 0.65(1) for the dodecagonal tiling (D = 6). Nev- 
ertheless, we emphasize that even for these codimensions, 
the value of the diffusion exponent is not constant over 
a wide time range and actually slightly decreases when 
the time increases. This leads us to conclude that either 
the asymptotic regime is not reached, or the long time 
behavior is not "strictly" power-law like. 
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For higher codimension, the situation is more compli- 
cated: no power-law regime can be clearly distinguished 
since j3 strongly depends on the time range and decreases 
when the time increases. 
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FIG. 3: Log-log plot of AR 2 {t)/t for high codimension 
quasiperiodic tilings. 

To determine whether it converges towards a asymtotic 
value (at large time) or not, we have focussed on the 
behavior of AR 2 (t)/t which gives us an indication on 
the super- or sub-diffusive character of the propagation. 
This is an important issue since in such two-dimensional 
quasiperiodic systems, no sub-diffusive propagation has 
ever been reported. As it can be seen in Fig. [3| although 
no power-law behavior can be extracted in the time range 
investigated, a strong indication of a diffusive asymptotic 
regime ( lim AR 2 (t)/t = Cte) is provided. 

t — >oa 

These results show that the codimension is an impor- 
tant parameter for the wavepacket propagation. Indeed, 
notin g th at for the two-dimensional generalized Rauzy 
tiling [lj] (D = 3) one obtains an exponent (3 = 0.95(1), 
it is clear that the higher the codimension, the lower /3. 
This can be easily understood in terms of the complexity 
of the structure for the following reasons. An impor- 
tant characteristic of quasiperiodic tilings is the so-called 
repetitivity property initially studied in two dimensions 
by Conway. In the class of tilings considered here, it 
states that for any local environment of linear typical 
size L, a similar environment can be found at a distance 
£ ~ L z with z > 1. In the simplest case D = 4, 5 one 
has z = 1. Moreover, when the codimension increases, 
the number of possible edge orientations at each node is 
lower or equal to 2D and for a given L, £(L) also grows. 
Consequently, the specific effects of the quasiperiodic or- 
der are expected to decrease when the codimension in- 
creases. Note that in periodic system, £ is a constant and 
in disordered systems, it grows exponentially with L. 

Nevertheless, these tilings are quasiperiodic and this 
special order is responsible, at least in low codimension, 
for a superdiffusive regime Q3 > 1/2). For higher codi- 
mension the quasiperiodicity becomes irrelevant and the 



exponent f3 seems to converge towards 1/2 which is the 
value expected for a weakly disordered system in two 
dimensions (for example a periodic lattice with random 
impurities or random magnetic fluxes [15j, or a disor- 
dered quasiperiodic Fibonacci lattice 16] ) where such a 
regime is expected when V AR 2 is lower than the localiza- 
tion length. Unfortunately, since the asymptotic regime 
is not reached with the system size studied here, we can- 
not estimate the codimension dependence of the diffusion 
constant. However, the results shown in Fig. [3] indicate 
that this constant does not seem to simply decrease when 
the codimension increases. 

If this analysis is correct, destroying quasiperiodicity 
in these tilings should also lead to a diffusive regime. 
To check this assumption and to corroborate our previ- 
ous statement, we have studied the quantum dynamics in 
random tilings obtained from the quasiperiodic tilings by 
making elementary phason flips (about 10 11 ) as shown in 
Fig. ^for D = 4. The random flip Markov chain has been 
studied into detail in the cases D — 3 and D = 4 (see 
[l7|). and the number of flips needed to get a random 
tiling uncorrelated to the initial quasiperiodic tiling is 
indeed of order N 2 . Furthermore, this number has been 
computed numerically in higher codimension tilings and 
is, in fact, always smaller than N 2 /2 [l^. In the present 
case, this leads to an upper bound of 5.10 11 flips. How- 
ever, we have checked for D = 4, that the spreading in 
tilings with 10 11 and 5.10 11 flips have essentially the same 
diffusion exponent. 
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FIG. 4: Log-log plot of AR 2 (t)/t for random tilings. 

As shown in Fig. 0] although no power-law regime 
can be clearly extracted, the propagation is very simi- 
lar to those displayed in Fig. |3l for high codimension 
quasiperiodic tilings. For the octagonal tiling (D = 4), 
this result refutes the conclusion of Passaro et al. claim- 
ing that the diffusion exponent increases when a phason 
disorder is introduced. This contradiction is due to the 
short time scale investigated in Ref. Indeed, we also 
observe that for t < 20, the propagation is faster in the 
D = A random tiling than in the D = 4 quasiperiodic 
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tiling. However, at larger times, the motion in the ran- 
dom tiling seems to become diffusive whereas, in the oc- 
tagonal tiling, it remains superdiffusive with ft = 0.77(1) 
[19( . We have also observed a similar crossover for higher 
codimcnsion (D = 5,6) (see Fig. |SJ. 
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FIG. 5: Log- log plot of /i = AR^ asi / AR^ adoln as a function 
of time. 

A possible explanation of these two regimes is that, 
at short times, the local quasiperiodic order is likely res- 
ponsible for coherent quantum interferences (absent in 
a random tiling) that slow down the dynamics. In this 
regime, the topological disorder thus favors, in a sense, 
the propagation. At large times, the influence of the 
quasiperiodicity remains the same but a weak localiza- 
tion process becomes dominant for the random tiling and 
leads to a diffusive regime. Of course, for larger codimcn- 
sion, the situation is different since the asymptotic regime 
is likely the same both in the quasiperiodic and in the 
random case. Thus, one expects /i = Ai?^ uaBi /Ai?^ ndom 
to become constant at large time and it is actually the 
case (see Fig. |S]for D = 10). 

This study raise interesting questions concerning the 
nature of the spectrum in quasiperiodic and random 
tilings. A simple way to characterize the spectral mea- 
sure of a state \tp) is to compute the averaged autocorre- 
lation function 

C(t) = - - f : dt'P(t')~t- Q , (4) 
1 Jo 

where P(t) = \(ifj(0)\ip(t))\ 2 . If a = 1, the spectral mea- 
sure is absolutely continuous, if a = 0, it is pure point, 
and if < a < 1 it is singular continuous. As previously, 
for a random phase state, a gives an average information 
over the total density of states. Unfortunately, in the 
tilings studied here, there exist strictly localized states 
known as confined states due to specific local environ- 



ment (see Ref. HilEl for the case D = 5). These states 
eventually leads to a — provided the initial random 
phase state has a nonzero overlap with them. The ana- 
lyzis of the spectral measure with respect to the codimen- 
sion thus requires to get rid of these confined states. Fur- 
ther investigations are needed to clarify this point which 
is certainly of great interest in the understanding of the 
electronic properties of quasiperiodic systems. 

We would like to thank M. Baakc and M. Duneau 
for fruitful remarks about the repetitivity property of 
quasiperiodic tilings. 
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